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D-dimensional constrained systems are studied with stochastic Lagrangian and 
Hamiltonian. It is shown that stochastic consistency conditions are second class constraints 
and Lagrange multiplier fields can be determined in (D+l)-dimensional canonical formula- 
tion. The Langevin equations for the constrained system are obtained as Hamilton's equa- 
tions of motion where conjugate momenta play a part of noise fields. 

The stochastic quantization scheme was first proposed by Parisi and Wu [11 as an al- 
ternative quantization method for gauge theories in 1981. After a few years Namiki et 
tt/.' 31 improved the method so as to handle constrained systems and recently more rigorous 
discussion was shown! 41 In their methods Lagrange multiplier fields can be determined under 
stochastic consistency conditions (the consistency conditions for the fictitious time). If we 
extend the stochastic quantization method to phase space, we can determine the Lagrange 
multipliers accompanied with the first class constraints by almost the same method.' 51 They 
seem strange because we cannot determine them in canonical formalism. In this letter we 
introduce the stochastic Lagrangian' 61 and Hamiltonian for constrained systems and treat 
them within a (D+l) -dimensional canonical framework. It is shown that stochastic consis- 
tency conditions are second class constraints and we obtain the Langevin equations for the 
constrained system as Hamilton's equations of motion. 



We consider a system described by the action S in D-dimensional Euclidean space-time. 
In the stochastic quantization scheme the Langevin equation is an important tool, which is 
written in the discretized form as 



5 S 

d qi = -—dt + dWi. (i = l,---,N) (1) 

oqi 

Here t is the fictitious time and dWi is the Wiener process whose expectation values are 
defined as 

{dWi(t)dWj(t') ■ ■ ■) = ^-J D(dW)dWi{t)dWj(t') ■ ■ -expj- J dxJ2 L ( dw ( T m))}, (2) 



where n is a normalization factor and L(dW) is stochastic Lagrangian. If the system has 
constraints 

F a = 0, (o=l,.--,M; N>M) (3) 
we define the stochastic Lagrangian as 

L{dW) = -^dWidWi + \ a dF a , (4) 

where A a is a Lagrange multiplier field. In the stochastic Lagrangian (4) we have changed 
the constraints (3) into their differentials. It is quite natural in this formulation because the 
order of the stochastic Lagrangian is dt and we rewrite it with the Langevin equations, i.e. 
stochastic differential equations, as shown in the following. 

When we handle stochastic differential equations, we should make choice of their calcu- 
lation rules, that is, adopt either Ito's calculation rule or Stratonovich's one in a general 
way. If we adopt the latter rule, we must bother about the Jacobian factor which appears 
when we change the integration variables in the generating functional. On the other hand 
the Jacobian factor is a unity if we calculate with Ito's rule) 61 therefore we will use it in the 
following calculation. In the Ito calculus we have Ito's formula 

(df) = (^-d qi ) + h^-I-dWidWj) (5) 
dqi 2 dqidqj 

where / is an arbitrary functional of qiS. 
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To change the integration variables we insert a unity 



/A Q 
D(q)5{d qi + —dt - dWi} 
oqi 

into RHS of (2) and integrate over dW . Then the distribution functional Z reads 

Z = j D(q)exp{- J rfx^L(g)}, (6) 



where 



L{q) = ^(d« + f^-dt) 2 + X a dF a , (7) 



d JjL dq . + I . d2Fa 
dqi 1 2 



dF a = -^d qi + - • -z-£rMij. (8) 



The second term in RHS of (8) is due to Ito's formula, My in which is a cit-order expectation 
value of dWidWj which is not 25{j because the Lagrange multiplier fields should contain a 
noise-like part. W We will compute it later. 

From the stochastic Lagrangian (7) we define stochastic conjugate momenta as 

dL I, dqi dL dF a 

» = Wi) = 2 ( Tit + + (9 ' a) 



The latter equation is a primary constraint, and hence the stochastic Hamiltonian H is 

if =Pidqi + n a d\ a - L + u a -K a dt 

=Pi(2pidt - ^-dt - 2X a ^^dt) - ^-(pidt - \J^dt)(pidt - X b ^-dt) 

5qi dqi dt dqi dqi (10) 

, ( dF a^ , SS , ox dF b ^ 1 d 2 F a „, -i 
- A {-7^(2pidt - — dt - Zh-g^dt) + - • ^-^-Mijj + u a ir a dt, 

where w a is the other Lagrange multiplier into which d\ a /dt is absorbed. 
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The stochastic differentials are defined using Poisson's bracket 



(AB\- dAdB dAOB OA OB OA OB 
dqi dpi dpi dqi d\ a dir a dn a d\ a 

as 

dA={A,H}. (fl) 

[3] 

The stochastic consistency condition of ir a pa thus reads 

d-K a = {lT a , H} 

= -g-(2 Pi dt - —dt) + - ■ Q-Q-Mij - 2——\ b dt (12) 

uq t uq t z, uq t uqj uq l uq t 

= (f) a paO. 

This is a secondary constraint. If we determine A& so as to satisfy it, we need no more 
constraints. We may identify <p a with dF a and n a with F a + C a where C a is an arbitrary 
functional independent of the fictitious time. Consequently we should choose the initial con- 
ditions as C a = to realize the constraints (3). The Poisson's bracket of the two constraints 
(9.b) and (12) is 

{n a ,4> b } = 2^^ = 2D ab . (13) 
dqi dq { 

Here we assume that D ab is not zero and has its inverse D ~^ . [3] Then the above equation 
shows that both of them are second class constraints. When we write the total Hamiltonian 
Ht as 

H T = H + v a 7r a dt + w a (j) a (14) 

where v a and w a are Lagrange multiplier fields, we can determine them with the help of the 
stochastic consistency conditions again. They read 

w a = 0, 
v a = l -D- ah l {h,H}. 



Accordingly the total Hamiltonian becomes 

1 



H T = H + ^D^{</> b ,H}n a . (15) 



4 



Next we derive Hamilton's equation of motion for qi. If we use Dirac's bracket 



[A, B] = {A, B} - {A, $ r }A-i{$ a , B} 



(16) 



with 



$ r = 7T a , a , 



A rs EE {$ r ,$ s }, 

we may replace the weak equality ~ by the strong equality =. In the following we always 
use Dirac's bracket instead of Poisson's bracket and rewrite the total Hamiltonian by the aid 
of the second class constraints (9.b) and (12) as 

Ht = PiPidt-pi—dt 
oqi 

-^{t (2wdt -^ t) + 5'^} (»> 

Using it we obtain Hamilton's equation of motion for qf. 

dqi = [qi,H T ] 

with 

which of course satisfies dF a = 0, and hence (N — M) independent equations exist. 

The equation (18) is Langevin equation if we may identify 2K{jpjdt with noise fields. 
In the following we show that is the case. In order to decompose the variables into con- 
straint variables and independent ones) 3141 we introduce a new set of variables {Q^}, (// = 
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1, 2, • • • , N), which are defined as 



(19) 



dqi % 

where 's are vielbein fields defined as 

^i = ^, e a ^D~ b ^, (a = N-M + l,...,N) (20.a) 

e1e a ,i = 0, (A=l,2,---,iV-M) (20.6) 

eA,i = (<T / S = e A A- (20.c) 

Here we have assumed that the metric g" 45 is non-singular. These definitions lead to following 
relations 

e ^i e v,i = Sfj, u , e ^i e fJ.,j — <%? (21. a) 

eVi = eie a ,i = e^e^k = e A ,ie a ,i = 0, (21.6) 
e i e A,j = K-iji e i e aj = -^y> (21- c ) 

Kije a 5 = Kijeaj = 0, (Sij - Kij)e A j = (Sij - Kij)e A j = 0, (21. d) 

det{e^) ^ 0. (21.e) 

The equations (21. e) shows that a manifold spanned by Q^s is equivalent to one by q^s. 

Assuming 2Kijpjdt to be the noise field and to give the same expectation value as dWi, 
which assumption will be confirmed later, we can derive the Langevin equations for the new 
variables. Taking account of Ito's formula (5), they are 

l<9e M - 

dQ» = e%dq i + --^M ij . (22) 

For the constraint variable Q a =F a , we have 

dQ a = (23.a) 
and the initial conditions should be chosen to be zero to realize the constraints (3),. On the 
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other hand, for the independent variables Q s, the Langevin equations are 

dQ A = -9 AB ^Bdt + e A (2 Pi dt) + ^e> y -. (23.6) 

Moreover if we decompose the vielbein e A as 

e^^EYi, (1 = 1,2,-. -N-M) (24.a) 

eVi = 5 /J , = E A E B j5 IJ = g AB , (24.6) 

where e I i only rotates the tangent space spanned by the independent variables, we may 
change the second term in RHS of the Langevin equation (23. b) into fields with N — M 
indices: 

e A (2 Pi dt) = E^dQ 1 , (25) 

where dVt 1 is defined as 

dtt 1 = e^^pidt), (26) 

which will be regarded as the Wiener process as shown in the following. 

We change the variables q^s in the distribution functional (6) into the new ones Q^s 

Z = J D(Q)det(e^)exp{- J dx^HQ)}, (27) 
and substitute zero for Q a, s. Then again we insert a unity 

/xq i f) p A 

D(dQ)det(E A )5{dQ A + g AB j^dt - E A dti - -^e^My} 

into the distribution functional (27) and integrate over Q A, s. Neglecting the normalization 
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factor, (27) becomes 

Z = J D(dtt)exp{- J dxJ2 L W}, (28) 



where 



L(dQ) = -^d^dVt!. (29) 



The expectation values are easily computed with the above distribution functional: 

(dQ 1 ) = 0, 

(drtdQ- 7 ) = 25 IJ dt, 

and 

My = 2(e A i e A ,j)dt. 
The Langevin equation (23. b) thus becomes 



xq f) p A 
dQ A = -9 AB j^dt + /:>/(>' + ^jfj<' ! ;dl. (30) 



which agrees with the Langevin equation for constrained systems obtained by the other 
method in the stochastic quantization scheme! 41 The last term is necessary for the general 
coordinate transformation covariance and, in perturbation theories, cancellation of some 
of divergent terms.' 41 The proper Fokker-Planck equation on (JV — M)-dimensional constraint 
surface is given by the ordinary prescription with the Langevin equations (23). It will be 
straightforward to extend the results to the phase-space stochastic quantization scheme. 

Finally we comment on the results when we obey Stratonovich's calculation rule. In this 
case the Langevin equation (30) should be a little modified.' 41101111 Nevertheless it, of course, 
does not mean the two rules lead different physical results. For example the Fokker-Planck 
equation, whose equilibrium distribution should be consistent with the path-integral form, 
is the same in both the cases. 
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